#include<iostream>
#include<algorithm>
#include<cmath>
using namespace std;
const int maxn=100;
struct point{
    double x,y;
    bool operator<(const point& p) const
    {
        return x<p.x;
    }
}p[maxn];
int n;
double dp[maxn][maxn];
int from[maxn][maxn];
double dis(const point&a,const point&b)
{
    return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}
// void print(int i,int j,int mark)
// {
//     // if(i<=1&&j<=2) return;
//     // if(i+2<=j) print(i,j-1,mark+1);
//     // else if(i+1==j) print(from[i][j],i,mark+1);
//     // else {
//     //     print(i-1,j,mark+1);
//     // }
//     // if(mark%2==0)
//     // cout<<from[i][j]<<" ";
// }
void print(int i,int j)
{
    if(i>=j)
    {
        if(i<=2&&i<=1) return;
        print(from[i][j],j);
        cout<<from[j][i]<<" ";
    }
    else{
        if(i<=i&&j<=2) return;
        print(i,from[i][j]);
        cout<<from[i][j]<<" ";
    }
}
int main()
{
    cin>>n;
    for(int i=1;i<=n;i++)
    cin>>p[i].x>>p[i].y;
    sort(p+1,p+1+n);
    dp[1][2]=dis(p[1],p[2]);
    from[1][2]=1;
    for(int i=1;i<=n;i++)
    {
        if(i>=2)
        dp[i][i+1]=dp[i-1][i]+dis(p[i+1],p[i-1]),from[i][i+1]=i-1;
        for(int k=1;k<=i-2;k++)
        {
            if(dp[i][i+1]>dp[k][i]+dis(p[i+1],p[k]))
            dp[i][i+1]=dp[k][i]+dis(p[i+1],p[k]),from[i][i+1]=k;
        }
        for(int j=i+2;j<=n;j++)
        {
            dp[i][j]=dp[i][j-1]+dis(p[j-1],p[j]);
            from[i][j]=j-1;
        }
        if(i==n)
        dp[i][i]=dp[i-1][i]+dis(p[i-1],p[i]),from[i][i]=i-1;
    }
    cout<<dp[n][n]<<endl;
    // print(n,n,0);
    print(n,n);
    // for(int i=1;i<=n;i++)
    // {
    //     for(int j=1;j<=n;j++)
    //     cout.width(10),cout<<from[i][j]<<" ";
    //     cout<<endl;
    // }
}